Superconducting phase transitions in ultrathin TiN films 
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We investigate transition to the superconducting state in the ultrathin (< 5nm thick) titanium 
nitride films on approach to superconductor-insulator transition. Building on the complete account 
of quantum contributions to conductivity, we demonstrate that the resistance of thin superconduct- 
ing films exhibits a non-monotonic temperature behaviour due to the competition between weak 
localization, electron-electron interaction, and superconducting fluctuations. We show that super- 
conducting fluctuations give rise to an appreciable decrease in the resistance even at temperatures 
well exceeding the superconducting transition temperature, Tc, with this decrease being dominated 
by the Maki- Thompson process. The transition to a global phase-coherent superconducting state oc- 
curs via the Berezinskii-Kosterlitz-Thouless (BKT) transition, which we observe both by power-law 
behaviour in current-voltage characteristics and by flux flow transport in the magnetic fleld. The 
ratio Tbkt/Tc foUows the universal Beasley-Mooij-Orlando relation. Our results call for revisiting 
the past data on superconducting transition in thin disordered films. 



It has long been known that in thin films the transition 
into a superconducting state occurs in two stages: first, 
with the decreasing temperature, the finite amplitude of 
the order parameter forms at the superconducting criti- 
cal temperature Tc, second, a global phase coherent state 
establishes at lower temperature Tbkt [the temperature 
of the Berezinskii-Kosterlitz-Thouless (BKT) transition] 
[l|, Q • Neither of these two temperatures results in sin- 
gularities in the temperature dependences of the linear 
resistivity R{T) and determination of Tc and therefore 
Tbkt from the experimental data is a non-trivial task. 
Fortunately, Tc can be determined by juxtaposing mea- 
sured R{T) with the results of the theory of supercon- 
ducting fluctuations (SF) in the region T > Tc, as it 
enters the theory as a parameter (see Ref. Q for a re- 
view). Tbkt can be determined either by the analysis of 
the power-law behaviour of current-voltage {TV) char- 
acteristics at zero magnetic field or by the change of the 
curvature of the magnetoresistance from the convex to 
the concave in the magnetic field perpendicular to the 
film plane 



Although superconductivity in ultrathin films has long 
been a subject of active investigations, only in the 
present work became it possible to comprehensively take 
into account all quantum contributions to conductivity 
(QCC). Our analysis is based on the recent theoretical 
advance [l^l offering a complete description of fluctua- 
tion superconductivity, which is valid at all temperatures 
above the superconducting transition. We demonstrate, 
in particular, that the omission of the Maki-Thompson 
contribution leads to the incorrect values of Tc- Our 
results thus call for revisiting the previously published 
data on the superconducting transition in thin disordered 
films. 

The data presented below are taken for thin (< 5nm) 



TiN films formed on a Si/Si02 substrate by atomic layer 
deposition and fully characterized by high-resolution 
electron beam, infra-red [l4| , and low-temperature scan- 
ning tunnelling spectroscopy [isl [l^. The films are 
identical to those that experienced the superconductor- 
insulator transition (SIT) (l7l - [l9t and have the diffusion 
constant D > 0.3 cm^ /sec and the superconducting co- 
herence length ^(i(0) 8nm. The samples were pat- 
terned into bridges 50 fim wide and 250 fim long. Trans- 
port measurements are carried out using low-frequency 
ac and dc techniques in a four probe configuration. 

We start with the fundamentals of the theory of quan- 
tum contributions to conductivity. A notion of quan- 
tum corrections to conductivity can be traced back to 
the work [20], where small corrections to the diffusion 
coefficient were discussed. In the years that followed a 
breakthrough in understanding the nature and role of 
quantum effects was achieved [21-2^], crowned by for- 
mulation of the concept of quantum contributions to con- 
ductivity in disordered metals [ii-i^. A theory of QCC 
(see for review 0, [2^) rests on two major phenomena: 
First, the diffusive motion of a single electron is accompa- 
nied by quantum interference of the electron waves [weak 
localization (WL)] impeding electron propagation. Sec- 
ond, disorder enhances the Coulomb interaction. The 
corresponding QCC comprise, in their turn, two compo- 
nents: the interaction between particles with close mo- 
menta [diffusion channel (ID)] and the interaction of the 
electrons with nearly opposite momenta (Cooper chan- 
nel). The latter contributions are referred to as super- 
conducting fluctuations (SF) and are commonly divided 
into four distinct types: (i) Aslamazov-Larkin term (AL) 
describing the flickering short circuiting of conductivity 
by the fluctuating Cooper pairs; (ii) Depression in the 
electronic density of states (DOS) due to drafting the 
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part of the normal electrons to form Cooper pairs above 
Tci (iii) The renormalization of the single-particle diffu- 
sion coefficient (DCR) [l^]; and last but not least (iv) the 
Maki- Thompson contribution (MT) [23, 24| coming from 
coherent scattering of electrons forming Cooper pairs on 
impurities. 

The total conductivity of the disordered system is thus 
the sum of all the above contributions added to the bare 
Drude conductivity Gq: 
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These contributions have their inherent temperature and 
magnetic field behaviours that govern the transport prop- 
erties of disordered systems and strongly depend on di- 
mensionality. We will focus on quasi-two-dimensional 
(quasi-2D) systems since our films fall into this cate- 
gory. In films the notation G in the Eq. ([ij refers to the 
conductance rather than to conductivity. The system 
is quasi-2D if its thickness, d, is larger than both, the 
Fermi wavelength, Xp, and the mean free path, /, but is 
less than the phase coherence length, l^p, (responsible for 
WL effect) and the thermal coherence length responsible 
for the electron-electron interaction in both the Cooper 
channel and the diffusion channel {It 
that is, Xp^ I < d < l^, It- 
corrections can be written as: 
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FIG. L (Color online) (a) Temperature dependences of super- 
conducting fluctuation contributions to conductivity [in the 
units of Goo = eV(27r^ft)] The curves for AL, DOS, and 
DCR processes are universal functions of reduced tempera- 
ture t = T/Tc, the MT contribution is presented for 5 = 0.01 
and 0.05. The black solid lines are the sum of all SF contribu- 
tions Eq. ([2]). The inset shows the same total SF contribution 
on logarithmic scale, (b) Resistance vs. reduced temperature 
(see details in the text), (c) The set of the curves Tmax/Tc 
vs. S for different coefficients A from Eq. ([3)|. The circles rep- 
resent the measured Tmax and Tc and 5 obtained by fitting 
the experimental data (see the Table I and discussion in the 
text). 



AG 



WL 



AG'^ =GooAln[fcsTT/ft] 



A = ap + Aid , 



(3) 



(4) 



where Gqo — e^/(27r^?i), a — 1 provided the spin-orbit 
scattering is neglected (r^ <^ Tso), a = —1/2 when scat- 
tering is relatively strong (Tcp ^ r^o), p is the exponent 
in the temperature dependence of the phase coherence 
time T^p cx T~p, and Ajd is a constant depending on the 
Coulomb screening and which in all cases remains of the 
order of unity (s^l • At low temperatures where electron- 
electron scattering dominates. 
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yielding p = 1, R\s the sheet resistance. At higher tem- 
peratures where the electron-phonon interaction becomes 
relevant, p = 2. This concurs with experimental observa- 
tions [siHiil where 1 < p < 2, with p = 1 at T < lOK. 

As far as superconducting fluctuations are concerned, 
the complete comprehensive formula, which includes all 
the SF contributions and is valid at all temperatures and 
magnetic fields above the superconducting transition line 
Bc2{T), was recently derived in Ref. [13|. We do not re- 



present the results of the calculations in the zero mag- 
netic field in Fig.HJ The AL, DOS, and DCR contribu- 
tions are universal functions of the reduced temperature 
t ~ T/Tc- The corresponding tabulated expressions are 
plotted in Fig.[T]i. The peculiarity of the MT contribu- 
tion is that it depends on the phase coherence time T^p, 
which enters through the parameter = -Kh/ {^ksTcT^). 
The latter can be expressed through the conventional 
pair-breaking parameter 



5 = 7rh/{8kBTT^) 



(6) 



as 7<p 



= tS. Note, that if t„ cx T ^, then S becomes 



temperature independent and by using the expression for 
Tip from Eq. ([5|) , Eq. ^ can be rewritten as 



'-R/{16h) \n[Trh/{e^R)] 
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produce here this somewhat cumbersome expression, but 



The MT contributions for the two values of the parameter 
S = 0.01 and S = 0.05 are presented in Fig.[T^. 

The following comments are in order. The calculations 
clearly demonstrate that at temperatures T > 2Tc the to- 
tal fluctuation-induced conductivity with the great pre- 
cision merely coincides with that given by the MT term. 
In other words, in this temperature range the sum of con- 
tributions AG^^ + AG^°^ + AG^^"=0. Moreover, the 
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MT process dominates the fluctuation superconductivity 
down to temperatures t — I < S. At lower temperatures 
the AL contribution starts to become larger. The fact 
that MT is a leading process in a wide range of tempera- 
tures has already been emphasized in very early works by 
Maki I23II and Thompson [2J| and subsequent theoretical 
works [2a, [S^] ■ On the experimental side, numerous ex- 
perimental studies that demonstrate that the magnetore- 
sistance at T > Tc is caused mainly by the suppression 
of the MT process 31 - i^l support this conclusion. 

Next, it is the common view that the temperature 
range (T — Tc)/Tc = t — l where superconducting fluctu- 
ations are relevant is defined by the so-called Ginzburg- 
Levanyuk parameter, t — 1 ^Gi. Note that in the two- 
dimensional case Gi= e^R/{23h) and as follows from ([7]) 
Gi« d. Therefore Gi defines only a narrow vicinity of Tc 
where the AL term becomes dominant. What concerns 
the total SF contribution, it remains noticeable and pos- 
itive even at temperatures well above Tc, as it is clearly 
demonstrated by the inset in Fig.[T^ presenting the total 
SF contribution on logarithmic scale. 

As the measurable quantity is the resistance, rather 
than the conductance, it is instructive to recast the cal- 
culated QCC into the temperature dependence of the re- 
sistance, R{T). To this end, we use Eq. ([7]) which offers a 
reasonable estimate for Go = in Eq. ([T]). In Fig.[T)D 
we show the temperature dependence of the resistance 
calculated as R{t) = [J:AG^'\t) + Go]"^ for 6 = 0.01 
yielding Gq ^ — 147^1. We choose the temperature refer- 
ence point at r = 20Tc assuming that all the contribu- 
tions are practically zero at this temperature. Summing 
up all the QCC we took the coefficient A = 2 in Eq. ([3]) . 
One sees that although the SF contributions alone would 
have resulted in the monotonic behaviour of the resis- 
tance (with dR/dT > 0), the contributions from WL 
and ID processes make R{T) become non-monotonic and 
exhibit a maximum at the some temperature, T^ax, of 
about of few Tc. 

In Fig.[T}: we plot the ratio T,nax/Tc as function of 5 
for three most common experimental situations where 
A = 2>,2 and 0.5, corresponding to three sets of pa- 
rameters (a, p, A/13) in Eqs. and (g]): (1,2,1), (1,1,1), 
and (-1/2,1,1). The immediate important conclusion to 
be drawn from these plots is that the maximum in the 
R{T) dependences for thin superconducting films is al- 
ways present, even for the smallest values of the coeffi- 
cient A. However, in relatively thick films this maximum 
can be disguised by the classical linear drop of the re- 
sistance with decreasing temperature. The next observa- 
tion is that for all the realistic parameters of the film, the 
larger the resistance R [i.e. the larger 5 Eq. ([7])], the closer 
Tmax moves to Tc, maintaining, nevertheless, that the 
ratio Tmax/Tc remains always larger than 2. It is note- 
worthy that the maximum lies in the domain where the 
SF are dominated by the Maki-Thompson contribution 
and that the maximum itself arises from the competi- 



TABLE I. Sample characteristics: -R300 is the resistance per 
square at T = 300 K; Rmax is the resistance at the maximum 
point achieved at temperature T = Tmax', Tc is the critical 
temperature determined from the quantum contribution fits 
together with the pair-breaking parameter 6 and the coeffi- 
cient A in Eq.Q; Tgj^j, is determined from the power- law 
behavior of TV characteristics and Tgj^j, is found from the 
flux flow resistance data. Note that the sample SOI is the 
same sample #1 in ref. [4l|], where the higher Tc = 0.6 K was 
found since only the AL contribution was taken into account. 
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tion between the WL-I-ID and MT processes. In general, 
Tmax/Tc VS. 5 curves relate the quantity Tmax which is 
the only characteristic point in the R{T) dependence with 
the transition temperature Tc and as such can serve as 
a set of calibrating curves for the express-determination 
of Tc, since A can be estimated from the analysis of the 
resistance behaviour at high temperatures. 

Now we turn to discussion of our experimental results. 
Figure[2^ presents the temperature dependences of the 
resistance per square for four TIN films with different 
room temperature resistances. In all samples the resis- 
tances first grow upon decreasing the temperature from 
room temperature down, then reach the maximum value, 
Rmax, at some temperature T^ax (see Table I), and, fi- 
nally, decreases, with Tmax being approximately three 
times larger than the temperature where R becomes im- 
measurably small. Before fltting the data with the theory 
of QCC, let us verify that the films in question are in- 
deed quasi-two-dimensional with respect to the effects of 
the electron-electron interaction. We find It ~ 2.5 nm at 
T = 300K {It > 12nm at T = 10 K). Therefore the con- 
dition of quasi- two-dimensionality, d < It, is satisfied at 
all temperatures down from 300 K. That is why the tem- 
perature behaviour of the conductance follows the loga- 
rithmic temperature dependence in accord with Eq. ([3]), 
see Fig.l^b. Solid lines in Fig.l^h-jb account for all quan- 
tum contributions. 

The fitting remarkably captures all major features 
of the observed dependences: their non-monotonic be- 
haviour, the position and the height of Rmax, and the 
graduate decrease in the resistance matching perfectly 
the experimental points down to values R <^ Rmax (with- 
out any additional assumptions about mesoscopic inho- 
mogeneities 42, i^]). We were using three fitting param- 
eters, S, A, and, Tc (the values providing the best fit are 
given in the Table 1). It is noteworthy that while varying 
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FIG. 2. (Color online) (a) Resistance per square vs. tempera- 
ture for four TiN film samples listed in Table I. Solid lines: fits 
accounting for all the corrections. Dashed lines (marked as 
WL+ID): separate contribution of the sum of weak localiza- 
tion and interaction in the diffusion channel to the resistance 
of the samples SOI and S15. Dotted lines (SF): contribu- 
tion of superconducting fluctuations, (b) The same data as 
in (a), but extended to room temperatures and re-plotted as 
the dimensionless conductance G/Goo- The semi-logarithmic 
scale representation reveals logarithmic decrease of the con- 
ductance with temperature owing to WL and ID effects, (c) 
Reduced conductance {R~^ — R~^^)~^ / Ral {Ra\ ~ e^/(16ft) 
and Rn = Rmax) vs. T. The linear fit to the AL expression 
(solid lines) are often used for the determination of Tc and 
generally gives incorrect (much too high) values of the criti- 
cal temperature: Td = 1.470 K, and 3.714 K for samples SOS 
and S04, respectively (marked by arrows). The correct values 
of Tc (presented in Table I) are in the rectangles and marked 
by the vertical bars. 



d and A significantly shifts the temperature position and 
the very value of Rmax, it does not change noticeably 
the position of Tc- It reflects the fact that AG^^ does 
not depend on the pair-breaking parameter S in the close 
vicinity of Tc (see inset to Fig.[T^ where the curves for 
different 6 merge). A cross-check of the validity of the 
extracted values of S and A is achieved by taking the val- 
ues of Tjnax/Tc and S found from fitting and superposing 
them on Fig.[T];. We see that the points fall between the 
lines corresponding to A = 2 and A = 3 in a nice accord 
with A ~ 2.6 obtained from the full description of R{T). 

Having completed the full description of the experi- 
mental data within the framework of a general theory 
of superconducting fluctuations 13], it is instructive to 
review the approaches for inferring Tc from the experi- 
mental data that were frequently used in the past. First, 
we find that Tc Hes at the foot of the R{T) curve where 



R{T) ~ (0.08^0.13)i?„aa;- Therefore, the determination 
of Tc as the temperature where R{T) drops to 0.5i?Ar 
(let alone to 0.9i?Ar) significantly overestimates Tc- An- 
other frequently used procedure [lo| is based on the as- 
sumption that the effect of quantum corrections can be 
reduced to the AL term only, i.e. that the resistance 
obeys the relation R^^ = R^^ + R']^\/{T/Tc - 1), where 
= e^/{im) = 1.52 • 10"^ f^"^ This implies that 
there would have existed the range of temperatures near 
Tc where the plot [{R'^ - R^^)'^ / Ral] vs. T could 
have been approximated with a straight line with the 
slope = 1. The intersection of this line with the T-axis 
would have defined Tc. Utilizing this approach we plotted 
in Fig.[2t the data for two of our samples, as an exam- 
ple, and indeed find such a linear dependence for each 
sample (shown by the solid lines) , yielding temperatures 
of the intersections marked as Tci- One sees, however, 
that this procedure gives much too high values for the 
superconducting critical temperatures. 

We now turn to the determination of Tbkt in our 
samples. There exists two distinct methods for finding 
Tbkt from transport measurements. The first utilizes 
power-law fits to the I-V characteristics of the film for 
which the switch from V I to V (x occurs at the 
transition f^. Figures[3^,b show the typical sets of the 
I-V curves for our samples at different temperatures near 
and below Tbkt as log-log plots which indeed represent 
V (X 1°' behaviour, with a rapidly growing in a narrow 
temperature window, characteristic to the BKT transi- 
tion (Fig.|3t). The values of T^^^j, determined as tem- 
perature where a = 3 are presented in the Table 1. Note 
that even at T > Tbkt while ^ cx / at low currents, 
the I-Vs become strongly nonlinear at elevated currents 
showing the characteristic rounding. In contrast to that 
at T < Tbkt there is an abrupt voltage jump termi- 
nating the power-law behaviour at a certain well defined 
current. We attribute this jump to a heating instability 
in the low-temperature BKT phase, the nature of which 
along with the formation of the critical current in the 
BKT state will be discussed in a forthcoming publica- 
tion 

The second technique for determini ng Tbk t involves 
the use of flux flow resistance data [3. [jl. [lol - [l3 | . Figure[3ji 
shows a family of magnetoresistance isotherms typical for 
all our samples. With increase of temperature the R{B) 
dependences progress from positive curvature through 
linear dependence to negative curvature. Below Tbkt 
field-induced free vortices not only contribute to the re- 
sistance due to their own motion, but screen antivortices 
helping to dissociate vortex-antivortex pairs. This results 
in a superlinear response to the applied field. Above 
Tbkt the vortex-antivortex pairs are unbound and all 
thermally induced vortices contribute to the resistance. 
Field-induced vortices annihilate with some fraction of 
antivortices, effectively reducing the number of vortices 
participating in the fiux-fiow, giving rise to a sublinear 
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FIG. 3. (Color online) Current-voltage characteristics on a 
double-logarithmic scale for samples S15 (a) and S04 (b) at 
different respective temperatures listed in the panels. Solid 
lines indicate the slopes corresponding to different values of 
power Q on the V oc 1°". (c) a{T) for sample S15 in a full- 
temperature scale showing that the development of the power- 
law behaviour occurs over a very narrow temperature interval. 
The inset presents the same data on a magnified temperature 
scale, emphasizing the jump from a = 1 to a = 3 at the tran- 
sition, (d) Resistance isotherms vs. magnetic field for sample 
S04 measured at temperatures 10 mK apart. The dashed lines 
are the tangents to lowest and highest isotherms stressing the 
evolution from superlinear to sublinear behaviour. 



response to the applied field. The temperature of the 
linear response indicates thus Tbkt- The correspond- 



are listed in Table L 



Note that while T^j^j, and T^^kt ^'"S very close to each 



ing temperatures denoted as Tj^^j, 

'I) 

other, the former appears to be slightly lower than the 
corresponding values of T^^j,. 

The summary of our results is presented in Fig. 21 
showing the ratio of Tbkt/Tc vs. i?3oo and Rmax- Irre- 
spectively to the choice of the resistance, the more disor- 
dered the film (i.e. the closer the film is to the SIT), the 
farther apart Tbkt from Tc is. Using the dirty-limit for- 
mula which relates the two-dimensional magnetic screen- 
ing length to the normal-state resistance Rn, Beasley, 
Mooij, and Orlando (BMO) [H (see also ^) proposed 
the universal expression for ratio Tbkt/Tc 
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FIG. 4. (Color online) Plot of Tbkt/Tc as a function of 
R* for TiN films. The triangles correspond to R* = i?3oo 
and squares stand for R* = Rmax giving better agreement 
with theory. Open symbols represent Tbkt determined from 
the power-law behaviour of I-V curves; closed symbols corre- 
spond to Tbkt obtained from fiux fiow resistance data. The 
lines represent BMO functions ^ for different values of /3 in 
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where A(T) is the temperature dependence of the super- 
conducting gap and parameter (3 = A(0)/(A;bTc). The 
ECS theory predicts value /3 = 1.76. The BMO function 
for this value of (3 is shown by the solid line in Fig. |4] and 
correctly describes the qualitative tendency of decreasing 
the Tbkt/Tc ratio with increasing disorder. The quan- 
titative comparison however encounters some problems. 
The source of one of them is that since disordered films 
exhibit strong temperature dependence of the resistance 
the choice of what value is to be taken as the normal-state 
resistance, Rn in Eq. ([8]), is not a priori obvious, see the 
extensive discussion of the ambiguity of the experimen- 
tal definition of Rn in [iHl. Indeed, we see that choosing 
i?30o as Rn, we obtain Tbkt/Tc dropping much faster 
that what is predicted by the BMO formula. At the same 
time taking Rmax as normal state resistance yields a re- 
sult more close to theory, see Fig. 21 Another source of 
discrepancy is that the BMO function contains the ra- 
tio A(0)/Tc. As it is shown in ref. [11], where the TiN 
films close by the parameters to those investigated in the 
present work were studied, this ratio is unusually large 
as compared to its BCS value and grows on the approach 
to the SIT. The dashed lines in Fig. H show the BMO 
function for elevated values of /3 and demonstrate that 
choosing values of /? corresponding to increasing disorder 
can improve the agreement between theory and experi- 
ment. 

In conclusion, we demonstrated that the temperature 
dependence of the resistance of quasi-two-dimensional 
superconducting films, including its non-monotonic be- 
haviour and the significant broadening of the transition 
is perfectly described by the theory of quantum contri- 
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butions to conductivity. The analysis based on careful 
account of all contributions enabled a precise determina- 
tion of the superconducting transition temperature. We 
found that the transition to the global phasc-cohcrcnt su- 
perconducting state occurs via the Berezinskii-Kosterlitz- 
Thouless transition, and that the ratio Tbkt /^c follows 
the universal Bcaslcy-Mooij-Orlando relation upon an 
appropriate choice of the normal state resistance and tak- 
ing into account the non-BCS A(0)/Tc ratio in disordered 
films. 
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